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HIGHLIGHTS 


•  The  proposed  system  protects  solar-collectors  from  overheating. 

•  The  system  combines  thermoelectric-self-cooling  technology  and  thermosyphon  effect. 

•  Heat  transfer  coefficients  reach  10.5  W/(m2  K),  75%  higher  than  a  static  dissipator. 

•  Proposed  measures  are  due  to  increase  this  value  up  to  15.23  W/(m2  K). 

•  No  electric  power  is  required  in  any  design. 
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Highly  promoted  by  the  European  Union  Climate  and  Energy  Package  for  2020,  solar  collectors  stand  out 
as  the  most  promising  alternative  to  meet  water  heating  demands.  One  of  the  most  limiting  problems  in 
these  systems  involves  the  overheating  of  the  working  fluid,  resulting  in  rapid  fluid  degradation,  scaling 
and  premature  component  failure. 

This  paper  presents  the  computational  design  of  a  zero-power-consumption  system  that  combines 
thermoelectric-self-cooling  technology  and  thermosyphon  effect  to  dissipate  the  excess  heat  from  a  real 
solar-collector  installation.  Thermoelectric  self-cooling  is  a  novel  thermoelectric  application  proven  to 
enhance  the  heat  dissipation  of  any  hot  spot  without  electricity  consumption. 

The  simplest  design  outperforms  currently-used  static  and  dynamic  dissipaters  for  overheating  pro¬ 
tection  in  solar  collectors,  since  it  increases  the  global  heat  transfer  coefficient  of  a  static  dissipater  by 
75%  and  requires  no  electricity.  Likewise,  the  final  design  presents  a  global  heat  transfer  coefficient  of 
15.23  W/(m2  K),  155%  higher  than  that  provided  by  static  dissipaters,  forming  a  reliable,  robust  and 
autonomous  system  that  stands  out  as  a  promising  alternative  to  prevent  the  overheating  of  solar 
collectors. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  energy  consumption  of  buildings  is  one  of  the  main  topics 
addressed  by  the  European  Union  Climate  and  Energy  Package  for 
2020.  This  is  no  surprise,  since  figures  indicate  that  this  sector  ac¬ 
counts  for  about  40%  of  the  primary  energy  consumed  in  the  Eu¬ 
ropean  Union  and  creates  over  40%  of  its  carbon  dioxide  emissions 
[1].  Deeper  into  these  figures,  a  quarter  of  this  energy  -around 
170  Mt  of  oil  equivalent  per  year-is  devoted  solely  to  meet  water 
heating  demands  [2].  The  Energy  Package  strongly  promotes  the 
use  of  renewable  sources  in  this  sector,  and  solar  energy  emerges  as 
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the  most  promising  alternative,  especially  in  southern  countries, 
where  the  high  levels  of  solar  radiation  allow  savings  of  about 
70-90%  of  total  water  heating  cost  [3]. 

In  a  solar  collector  system  [4],  solar  radiation  is  absorbed  and 
the  resulting  heat  is  used  to  warm  fresh  water  in  a  tank.  The  heat  is 
transferred  from  the  collectors  to  the  tank  by  a  fluid  (usually  a 
mixture  of  water  and  glycol)  that  circulates  by  natural  convection 
(passive  solar  systems)  or  forced  by  a  pump  (active  solar  systems). 
These  systems  are  designed  to  cover  almost  all  the  water  heating 
demands  of  a  building  in  summer,  and  from  40  to  80%  in  winter, 
providing  hot  water  at  around  60  °C  [3,4]. 

One  of  the  most  important  problems  of  these  systems  is  the 
overheating  of  the  working  fluid,  which  takes  place  when  this  fluid 
surpasses  its  critical  high-limit  temperature  and  starts  boiling. 
Repeated  overheating  causes  glycol  to  degrade  and  become  acidic, 
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Nomenclature 

Pc 

pressure  loss  in  the  collectors,  Pa 

Pp 

pressure  loss  in  the  connecting  pipes,  Pa 

A 

area,  m2 

PjSC 

pressure  loss  in  the  TSC,  Pa 

Ap 

aperture  area  of  the  solar  collector,  m2 

Pr 

Prandtl  number 

Cll 

main  factor  for  thermal  loss  of  the  solar  collector,  W / 

Peltier 

Peltier  effect 

(m2  K) 

Q 

heat  flow  rate  absorbed  by  the  fluid  in  the  solar 

a2 

secondary  factor  for  thermal  loss  of  the  solar  collector, 

collectors,  W 

W/(m2  K2) 

Re 

Reynolds  number 

C 

thermal  capacity,  J/K 

Rc 

thermal  resistance  of  the  cold  side  heat  exchanger, 

c 

specific  heat  of  the  fluid,  J(/kg  K) 

K/W 

cold 

cold  side  of  the  thermoelectric  pairs 

Rdis 

thermal  resistance  of  each  dissipater,  K/W 

DP 

diameter  of  the  connecting  pipes,  m 

Rh 

thermal  resistance  of  the  hot  side  heat  exchanger,  K/W 

Disc 

diameter  of  the  channels  of  the  TSC,  m 

Rc 

electrical  resistance  of  a  thermoelectric  module,  Q 

Du 

number  of  design  units 

S 

cross  section  of  the  channels,  mm2 

E 

electromotive  force  produced  by  a  thermoelectric 

Te 

environment  temperature,  °C 

module,  V 

Th 

average  fluid  temperature  in  the  channels  of  the  TSC,  °C 

f 

friction  factor 

Ti 

fluid  temperature  at  the  inlet  of  the  solar  collectors,  °C 

g 

gravity  of  earth,  m/s2 

T0 

fluid  temperature  at  the  outlet  of  the  solar  collectors,  °C 

H 

height  of  the  solar  collector,  m 

r 

temperature  in  the  time  instant  t  +  At,  °C 

h 

vertical  distance  between  the  collectors  and  the  TSC,  m 

ThomsonThomson  effect 

hh 

convective  coefficient  in  the  channels  of  the  TSC,  W/ 

t 

time,  s 

(m2  I<) 

U 

heat  transfer  coefficient,  W/(m2  I<) 

hot 

hot  side  of  the  thermoelectric  pairs 

V 

volumetric  flow  rate  of  the  fluid,  m3/s 

I 

solar  radiation,  W/m2 

z 

vertical  height  coordinate,  m 

J 

electric  current  through  the  thermoelectric  modules,  A 

AV 

electric  voltage  supplied  by  four  thermoelectric 

Joule 

Joule  effect 

modules  to  a  fan,  V 

I< 

coefficient  of  minor  pressure  losses 

At 

time  step  of  the  computational  model,  s 

k 

thermal  conductivity  of  the  fluid,  W/(mI<) 

a 

Seebeck  coefficient,  V/K 

LP 

length  of  the  connecting  pipes,  m 

P 

tilt  angle  of  the  solar  collector,  ° 

Ljsc 

length  of  the  channels  of  the  TSC,  m 

V 

thermal  efficiency  of  the  solar  collector 

M 

number  of  pairs  in  a  thermoelectric  module 

Vo 

optical  efficiency  of  the  solar  collector 

N 

number  of  thermoelectric  modules  installed  in  the  TSC 

V 

kinematic  viscosity  of  the  fluid,  m2/s 

n 

number  of  channels  of  the  TSC 

P 

density  of  the  fluid,  kg/m3 

P+ 

pressure  increase  provided  by  thermosyphon  effect,  Pa 

T 

Thomson  coefficient,  V/K 

P~ 

pressure  loss  in  the  closed  loop,  Pa 

accelerates  scaling  and  causes  premature  component  failure  [4,5]. 
The  overheating  of  the  working  fluid  can  take  place  under  electric 
power  failures,  circulator  pump  failures  and  maintenance  shut¬ 
downs.  More  often,  it  occurs  under  excess  solar  radiation,  espe¬ 
cially  in  summer,  when  the  system  collects  high  amounts  of  heat 
but  the  hot  water  demand  is  at  its  lowest  level. 

Nowadays,  there  is  a  wide  variety  of  methods  to  provide  over¬ 
heating  protection.  The  optimal  system  must  work  under  all  situ¬ 
ations,  including  those  where  no  electric  power  from  the  grid  is 
available.  Furthermore,  it  would  seem  quite  contradictory  from  an 
energy  efficiency  point  of  view  that,  under  overheating  conditions, 
not  only  the  energy  harvested  by  the  collectors  was  wasted,  but 
also  additional  electric  power  was  required  to  remove  this  excess. 

In  the  first  place,  the  most  novel  methods  involve  the  use  of 
scattering  layers,  thermotropic  layers  or  substances  with 
temperature-based  solubility,  to  regulate  the  absorbed  heat  in  or¬ 
der  to  meet  instant  hot  water  demands,  without  electric  power 
consumption.  Still  under  deep  investigation,  these  technologies  are 
due  to  cope  with  overheating  problems  in  the  future,  as  soon  as 
their  long-term  performance  is  proven  and  their  cost  reduced  [5]. 

In  the  second  place,  drain-back  systems  remove  the  fluid  from 
the  pipes  before  it  surpasses  its  critical  high-limit  temperature  [4]. 
In  the  non-electricity-consuming  version,  the  increasing  pressure 
in  the  pipes  forces  the  fluid  to  move  to  a  reservoir  tank.  Then,  under 
high  solar  radiation,  the  collectors  heat  up  quickly  and  reach  their 
equilibrium  temperature  (called  stagnation  temperature).  The 


main  problems  of  these  systems  are  caused,  firstly,  by  the  presence 
of  remaining  fluid  in  the  pipes  that  overheats  and  degrades,  and, 
secondly,  by  the  accelerated  ageing  of  the  mechanical  components 
under  repeated  stagnation  conditions. 

Finally,  many  solar  collectors  include  static  or  dynamic  dissi- 
paters  to  transfer  the  excess  heat  to  the  ambient.  In  the  non¬ 
electricity-consuming  version,  the  closed  loop  composed  of  the 
collectors,  the  static  dissipaters  and  the  connecting  pipes  is  iso¬ 
lated,  and  the  thermosyphon  effect  creates  the  force  required  to 
make  the  fluid  circulate  [6].  Flowever,  static  dissipaters  work  by 
natural  convection,  and  this  mechanism  limits  their  performance  to 
a  great  extent,  imposing  low  heat  transfer  coefficients  and  high 
sensibility  to  environment  conditions  [7]. 

A  previous  paper  [8]  shows  a  promising  new  alternative  for 
heat  dissipation,  called  thermoelectric  self-cooling  (TSC),  which 
joins  the  advantages  of  dynamic  and  static  dissipaters.  This 
application  uses  thermoelectric  generation  technology  to  improve 
the  heat  transfer  of  a  passive  cooling  system,  such  as  a  static 
dissipater,  without  electricity  consumption.  In  a  TSC,  several 
thermoelectric  modules  are  installed  between  a  hot  spot  and  a 
passive  cooling  system.  These  modules  transform  part  of  the  heat 
into  electricity,  which  is  supplied  directly  to  the  cooling  system.  As 
a  result,  the  cooling  system  becomes  active  and  its  heat  transfer 
coefficient  increases.  The  use  of  TSC  for  overheating  protection  has 
been  explored  in  a  previous  paper  [9],  presenting  a  dissipation 
system  that  increases  by  more  than  50%  the  global  heat  transfer 
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coefficient  provided  by  static  dissipaters  currently-used  in  solar 
collectors,  and  eliminates  the  influence  of  the  atmospheric  wind 
and  the  limitations  connected  to  natural  convection.  However,  this 
TSC  requires  that  the  circulation  pump  works  continuously,  which 
renders  it  useless  under  the  electric  power  failures  and  pump 
malfunctions. 

The  present  paper  proposes  the  idea  of  combining  TSC  tech¬ 
nology  and  thermosyphon  effect  to  provide  overheating  protec¬ 
tion  to  solar  collectors  for  water  heating.  The  main  objective  is  to 
develop  and  a  computational  design  of  a  reliable  and  autono¬ 
mous  zero-power-consumption  TSC  system,  intended  to  dissi¬ 
pate  the  excess  heat  in  a  real  solar-collector  installation,  and 
assess  its  performance.  To  this  end,  the  paper  is  set  as  follows: 
Section  2  describes  the  real  solar-collector  installation  and  de¬ 
fines  the  overheating  conditions;  then,  Section  3  introduces  the 
main  components  of  the  TSC  and  describes  its  performance;  after 
that,  Section  4  explains  the  methodology  used  to  conduct  the 
computational  design  and  introduces  the  expressions  that  define 
the  thermosyphon  effect  and  the  simulation  of  the  TSC.  Finally, 
the  results  and  corresponding  analysis  are  presented  in  Section 
5,  whereas  Section  6  summarizes  the  main  conclusions  of  the 
work. 

2.  Solar  collector  installation 

The  solar-collector  system  is  installed  in  a  one-family  house  in 
the  north  of  Spain.  It  is  composed  of  two  flat-plate  collectors 
Astersa  AS  -  2.6  M,  set  in  parallel,  both  with  55°  of  tilt  angle,  and  a 
storage  tank  for  a  maximum  of  260  L  of  hot  water  at  60  °C.  The 
characteristics  of  this  solar  collector  are  presented  in  Table  1  [10]. 
The  thermal  fluid  is  a  mixture  of  water  (60%)  and  glycol  (40%)  that 
circulates  under  the  action  of  a  pump  with  40  W  of  electric  power 
consumption.  Eqs.  (1)— (4)  provide  respectively  the  temperature- 
dependant  density,  specific  heat,  thermal  conductivity  and  kine¬ 
matic  viscosity  of  the  fluid  11  ]. 


1055  —  0.70T 

(1) 

3612  +  3.75T 

(2) 

(350  +  1.80T) 

10 3 

(3) 

(3.02  -  0.03T 

)l(T6 

(4) 

The  most  unfavourable  overheating  conditions  occur  when  the 
storage  tank  is  full  with  water  at  60  °C  and  there  is  no  demand,  so 
that  all  the  heat  harvested  by  the  collectors  must  be  evacuated  by 
the  dissipation  system.  Moreover,  the  maximum  solar  radiation  in 
the  zone  (1000  W/m2)  takes  place  and  combines  with  the 
maximum  environment  temperature  (35  °C)  and  negligible  wind 
speed  12].  To  prevent  the  overheating  of  the  thermal  fluid,  the 
system  is  endowed  with  two  fin-tube  static  dissipaters,  set  in  par¬ 
allel  and  installed  after  the  collectors.  Each  dissipater  presents  6 
finned  tubes  with  length  1000  mm  and  diameter  16  mm,  outer  heat 
transfer  area  of  4.73  m2  and  global  heat  transfer  coefficient  (under 


Table  1 

Characteristics  of  the  solar  collector  Astersa  AS  —  2.6  M. 


Aperture  area  Height  Width  Number  of  channels 

2.54  m2  2.18  m  1.259  m  10 

Optical  Main  factor  for  Secondary  factor  for  Diameter  of  the 

efficiency  thermal  loss  thermal  loss  channels 

0.748  3.44  W/(m2 1<)  0.015  W/(m2  K2)  8  mm 


no-wind  conditions)  of  5.96  W/(m2 1<)  [13].  The  horizontal  area  and 
volume  occupied  by  the  two  dissipaters  reaches  0.48  m2  and 
0.019  m3  respectively.  When  the  fluid  leaving  the  collectors  reaches 
90  °C,  it  is  diverted  from  the  original  circuit  and  forced  to  flow 
though  the  dissipaters,  losing  the  heat  previously  absorbed.  Then,  it 
returns  to  the  original  circuit,  crosses  the  tank  and  goes  back  to  the 
collectors.  The  circulation  pump  must  work  continuously  during  all 
this  process. 

3.  Thermoelectric  self-cooling  system 

Fig.  1  presents  the  proposed  TSC  for  overheating  protection, 
which  comprises  standard  thermoelectric  modules  Kryotherm 
TGM-287-1.0-1.5  [14]  connected  electrically  in  series  and  thermally 
in  parallel,  that  harvest  heat  from  the  fluid  and  transform  it  into 
electricity.  The  hot  side  heat  exchanger  (HSHE)  connects  the  fluid 
and  the  hot  side  of  the  modules,  and  is  composed  of  several 
aluminium  plates  with  circular  channels.  Finally,  the  cold  side  heat 
exchanger  (CSHE)  transfers  heat  from  the  cold  side  of  the  modules 
to  the  environment,  comprising  several  finned  dissipaters  and  axial 
fans  energized  by  the  modules. 

When  the  fluid  temperature  at  the  outlet  of  the  collectors  rea¬ 
ches  90  °C,  it  is  forced  to  flow  through  the  channels  of  the  HSHE. 
This  heat  exchanger  warms  up,  inducing  a  temperature  difference 
between  the  ends  of  the  thermoelectric  modules.  Then,  the  mod¬ 
ules  provide  electric  power  to  the  fans,  forcing  air  to  flow  over  the 
dissipaters.  As  a  result,  the  thermal  resistance  between  the  fluid 
and  the  environment  decreases,  enhancing  the  global  heat  transfer 
coefficient.  The  fluid  cools  down  as  it  crosses  the  TSC  and  returns  to 
the  collectors. 


0.15  m 


0.10  m 


DISSIPATOR 

MODULES 

THERMAL 

FLUID 

METALLIC 

PLATES 

FAN 


FLUID 

INFLOW 


Fig.  1.  Front  and  side  view  of  a  design  unit  of  the  thermoelectric  self-cooling  system. 
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4.  Methodology 

4.1.  Thermosyphon  effect 

As  can  be  seen  in  Fig.  2,  once  the  fluid  is  enclosed  in  the  dissi¬ 
pation  loop,  the  thermosyphon  effect  establishes  volumetric  flow 
rate  that  circulates  anticlockwise  [6].  This  volumetric  flow  rate  is 
such  that  the  pressure  increase  provided  by  the  thermosyphon 
effect  equals  the  pressure  loss  in  the  collectors,  channels  of  the  TSC 
and  connecting  pipes,  as  Eq.  (5)  indicates.  Both  terms  of  Eq.  (5)  are 
explained  in  subsections  4.1.1  and  4.1.2. 

p+(v)=p-(v)  (5) 

4.1.1.  Pressure  increased  by  thermosyphon  effect 

The  pressure  increase  in  the  fluid  is  provided  by  Eq.  (6).  All  the 
channels  of  the  TSC  are  horizontal  and  no  temperature  gradient  has 
been  considered  in  the  connecting  pipes. 

Hsen/3 

P+(v^)  =  f  {Pg)dz  =  (Hsen(3  +  h)gPli-  j  Pgdz-hgPTo  (6) 

0 

All  the  terms  are  already  fixed  in  previous  sections,  except  the 
fluid  temperature  at  the  collector  inlet.  Since  the  outlet  tempera¬ 
ture  is  set  to  90  °C  -which  is  the  highest  allowed  to  prevent  the 
overheating-the  inlet  temperature  is  calculated  with  the  thermal 
balance  in  the  collectors,  provided  by  Eq.  (7).  This  inlet  temperature 
and  also  the  final  value  of  the  pressure  increase  are  functions  of  the 
volumetric  flow  rate. 


Q  =  Vpc(T0-Ti)=AIv=AI 


VO 


/ 


a2((ri  +  r0)/2-re) 


i 


(7) 


Eq.  (6)  indicates  that  P+  increases  as  h  increases,  since  the  fluid 
density  at  the  inlet  temperature  is  always  higher  than  that  at  the 


outlet  temperature,  according  to  Eq.  (1),  so  the  vertical  distance 
between  the  collector  outlet  and  the  dissipation  system  can  be 
increased  if  needed.  However,  h  higher  than  zero  means  that  the 
dissipation  system  will  be  under  direct  solar  radiation,  reaching 
high  temperatures  while  being  in  standby.  If  so,  the  first  ounces  of 
fluid  might  boil  inside  the  dissipation  system,  causing  steam 
hammering  and  extremely  rapid  temperature  change  that  can  lead 
to  premature  ageing  and  failure.  To  prevent  this,  the  dissipation 
system  must  be  installed  either  behind  the  collectors,  thus  reducing 
P+  and  limiting  its  applicability  to  collectors  with  high  tilt  angle,  or 
at  higher  level  but  covered  with  a  protective  layer.  It  should  be 
noted  that  the  latter  would  be  useless  if  a  static  dissipater  was  used 
as  dissipation  system,  since  the  cover  would  be  an  obstacle  for  the 
air  stream  created  by  natural  convection,  reducing  even  more  the 
performance  of  the  static  dissipater  7].  However,  this  factor  has  no 
significance  if  a  TSC  is  used,  since  forced  convection  is  imposed  at 
the  outer  surface  and  the  air  stream  is  easily  oriented. 


4.1.2.  Pressure  loss  in  the  circuit 

Eq.  (8)  provides  the  pressure  loss  in  the  closed  loop,  which  is 
composed  of  that  in  the  collectors,  connecting  pipes  and  dissipa¬ 
tion  system. 

P  =^c  +  ^TSC  (8) 

First  of  all,  the  pressure  loss  in  a  collector  is  provided  by  the 
manufacturer  as  a  function  of  the  volumetric  flow  rate,  and  is 
presented  in  Table  2. 

Secondly,  Eq.  (9)  presents  the  general  expression  for  the  pres¬ 
sure  loss  in  a  circular  pipe.  It  is  clear  that  the  fluid  regime  in  the 
connecting  pipes  must  be  laminar;  otherwise,  the  pressure  loss 
would  be  so  high  that  the  thermosyphon  effect  might  not  suffice  to 
move  the  fluid.  Pipe  length  is  set  to  4  m,  which  seems  enough  to 
connect  the  collectors  and  the  dissipation  system,  whereas  pipe 
diameter  is  set  to  26  mm  to  ensure  laminar  flow.  Then,  friction 
factor  can  be  obtained  with  Eq.  (10)  as  a  function  of  the  Reynolds 
number,  provided  by  Eq.  (11).  To  conclude,  a  quite  conservative 
value  of  10  is  set  for  the  coefficient  of  minor  pressure  losses  in  the 
pipes,  associated  with  bends,  fittings,  valves,  etc  [6]. 


pf-(ftrp+K) 

\v2 

7T2  Dj 

(9) 

/  =  64 /Re 

(10) 

II 

4 

(ID 

Finally,  the  last  term  in  Eq.  (8)  represents  the  pressure  loss  in  the 
channels  of  the  TSC.  In  the  abstract,  the  volumetric  flow  rate  is 
different  in  each  channel.  However,  given  that  laminar  flow  must 
be  imposed,  this  inequality  will  be  insignificant  and  it  will  a  valid 
approximation  to  consider  that  equal  volumetric  flow  rate  circu¬ 
lates  through  all  the  channels  [6].  Then,  Eqs.  (9)— (11)  are  valid  to 
calculate  the  pressure  loss  in  a  single  channel,  wherein  the  length 
and  the  diameter  are  design  parameters,  the  volumetric  flow  rate  is 
the  total  volumetric  flow  rate  divided  by  the  number  of  channels, 
and  the  coefficient  of  minor  pressure  losses  are  set  at  2,  comprising 
two  flanged  tees  [6]. 


Table  2 

Pressure  loss  in  the  solar  collector  Astersa  AS  —  2.6  M. 

V  x  105  (m3/s)  1.11  1.39  2.22  3.33 

Pc  (Pa)  65  80  100  200 
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4.2.  Simulation  of  the  TSC 

A  computational  model,  created  specifically  to  simulate  TSC 
applications  15],  is  used  to  determine  the  system  performance. 
Based  on  the  finite  differences  method,  the  model  is  able  to 
simulate  the  thermoelectric  system,  including  the  influence  of 
thermal  and  electric  contacts  and  all  thermoelectric  effects  (See- 
beck,  Peltier,  Joule  and  Thomson)  with  temperature-dependant 
properties. 

The  computational  model  transforms  the  system  into  a  grid  of 
connected  nodes,  and  solves  the  non-linear  set  of  equations 
composed  of  the  thermoelectric  effects  and  the  three-dimensional 
Fourier's  law  in  implicit-finite-differences.  Eq.  (12)  provides  this 
expression  applied  to  a  node  “a”  connected  to  x  additional  nodes, 
whereas  Eq.  (13)  provides  the  thermal  resistances  between  this 
node  and  the  x  additional  ones,  where  A  stands  for  the  cross  area 
between  nodes  and  U  stands  for  the  corresponding  convective, 
conductive  or  contact  heat  transfer  coefficient.  As  can  be  seen,  the 
model  uses  the  temperatures  of  all  the  nodes  (T),  heat  flow 
rates (Q),  thermal  resistances  between  nodes  (R)  and  thermal  ca¬ 
pacities  (C)  to  calculate  temperatures  in  the  subsequent  time 
instant  (T) 

_l_(r'_7;)  +  ...+_L(7'_T;)+Qa=g(r'-rfl)  i= 

(12) 

Rai  =  X/uaiAa  1  =  (13) 

Inside  a  thermoelectric  module,  heat  is  absorbed  by  Peltier  ef¬ 
fect  at  the  hot  ends  of  the  M  pairs,  whereas  heat  is  generated  at  the 
cold  ends.  Eqs.  (14)  and  (15)  present  these  effects.  Likewise,  Eq.  (16) 
provides  the  heat  generated  by  Joule  and  Thomson  effects  along  the 
semiconductor  legs,  being  Ro  the  electrical  resistance  of  the  ther¬ 
moelectric  module  and  r  the  Thomson  coefficient,  derived  from 
Seebeck  coefficient  by  using  Eq.  (17).  Finally,  Eq.  (18)  provides  the 
electromotive  force  produced  by  a  thermoelectric  module. 

Q.Peltier,hot  =  —M/7" hotshot  (14) 


Eq.  (7)  provides  the  fluid  temperature  after  the  TSC  and  Eq.  (19) 
calculates  the  average  logarithmic  temperature  [7].  This  average 
logarithmic  temperature  is  used  as  heat  source  temperature  in  the 
simulations  with  the  computational  model. 


Th  —  Te  + 


Ti  -  To 


ln((Tf  -  Te)/(T0  -  Te)) 


(19) 


4.2.2.  Thermoelectric  modules 

Kryotherm  TGM-287-1.0-1.5  thermoelectric  module  comprises 
287  pairs  of  n-doped  and  p-doped  bismuth-telluride  legs.  Each  leg 
presents  base  area  of  1  x  1  mm2  and  length  of  1.5  mm.  They  are 
sandwiched  between  two  similar  plates,  each  of  them  with  base 
area  of  40  x  40  mm2  and  length  of  1.5  mm,  that  provide  electrical 
insulation.  The  number  of  thermoelectric  modules  installed  in  the 
TSC  is  a  design  parameter. 


4.2.3.  Thermal  resistance  of  the  hot  side  heat  exchanger 

The  thermal  resistance  of  the  HSHE  includes  three  effects:  heat 
convection  between  the  fluid  and  the  channel  wall,  heat  conduc¬ 
tion  through  the  aluminium  plates  and  thermal  contact  between 
the  plate  surface  and  the  modules.  Firstly,  this  last  term  has  been 
deeply  studied  and  common  values  in  the  literature  indicate  a  heat 
transfer  coefficient  of  about  21,000  W/(m2  I<),  corresponding  to  a 
thermal  resistance  of  0.03  K/W  for  a  thermoelectric  module  with 
40  x  40  mm2  of  base  area  16].  Secondly,  the  conductive  thermal 
resistance  can  be  dismissed,  given  the  high  thermal  conductivity  of 
aluminium  and  the  minimal  thickness  that  the  plates  should  have. 
Finally,  Eq.  (20)  presents  Sieder  and  Tate's  expression  to  calculate 
the  corresponding  convective  heat  transfer  coefficient  between  the 
laminar-flow  fluid  and  the  wall  of  a  short  channel  [7].  Reynolds  and 
Prandtl  numbers  are  respectively  provided  by  Eqs.  (11)  and  (21). 
The  number  of  channels  and  the  diameter  of  a  channel  are  also 
design  parameters.  Then,  Eq.  (22)  provides  the  thermal  resistance 
of  the  HSHE. 


hhDjsc 

k 


1.86 


(20) 


Qpeltier,cold  —  ^JTco\daco\d 


(15) 


(21) 


Qjoule, Thomson  —  JR  0  Mr7(Tj10t  Tco\^) 


(16)  ^ ^/hfj (htcDtsc^tsc) 


(22) 


r  =  Tda/dT  (17) 

E  =  M((a\ ^otThot  —  acold^cold)  —  T(^hot  —  ^cold))  0&) 

Fig.  3  shows  the  flow  chart  of  the  computational  model.  The 
inputs  are  the  temperature  of  the  heat  source  and  the  environment, 
the  thermal  resistances  of  the  heat  exchangers  and  the  number  and 
type  of  thermoelectric  modules,  selected  from  a  data  base  that 
provides  the  dimensional,  electrical  and  thermal  parameters.  The 
model  outputs  are  voltage,  electric  current  and  electric  power 
generated  by  the  modules,  efficiency  of  the  thermoelectric  con¬ 
version,  temperatures  and  heat  flow  rates  transferred  between 
components  and  generated/absorbed  by  Joule,  Thomson  and  Peltier 
effects.  Reported  deviations  between  experimental  and  simulated 
outputs  lay  between  ±12%  [15]. 

4.2.1.  Heat  source  and  environment 

As  Section  2  indicates,  the  ambient  temperature  is  set  at  35  °C, 
whereas  the  fluid  temperature  before  the  TSC  is  set  at  90  °C.  Then, 


4.2.4.  Thermal  resistance  of  the  cold  side  heat  exchanger 

The  thermal  resistance  of  the  CSHE  includes  the  thermal  contact 
between  the  modules  and  the  dissipaters,  the  heat  conduction 
through  the  dissipaters  and  the  heat  convection  between  the  dis¬ 
sipaters  and  the  environment.  This  last  term  is  not  obvious,  since  it 
depends  on  the  air  stream  forced  over  the  fins,  which  in  turn  de¬ 
pends  on  the  electric  power  supplied  by  the  modules  to  the  fans. 
Therefore,  the  thermal  resistance  of  the  CSHE  must  be  introduced 
as  a  function  of  this  electric  power. 

The  calculation  of  the  mathematical  function  that  relates  these 
two  variables  requires  the  use  of  computational-fluid-dynamics 
software  and  experimental  procedures  [15].  This  calculation  has 
been  carried  out  for  the  assembly  composed  of  four  thermoelectric 
modules  Kryotherm  TGM-287-1.0-1.5,  connected  electrically  in  se¬ 
ries  and  thermally  in  parallel,  attached  to  an  aluminium  finned 
dissipater  comprising  a  square  base  plate,  with  side  length  150  mm 
and  height  12  mm,  and  23  fins  with  dimensions  150  x  25  x  1.5  mm3. 
A  fan  Sunon  KDE1208PTS1  is  energized  by  the  four  modules, 
providing  forced  conversion  over  the  dissipater.  Eq.  (23)  shows  the 
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Fig.  3.  Flow  chart  of  the  computational  model. 


obtained  relation  between  the  thermal  resistance  of  this  dissipater 
and  the  voltage  supplied  by  the  four  modules  to  the  fan  [15].  Since 
the  CSHE  includes  several  of  these  assemblies,  the  thermal  resis¬ 
tance  of  this  heat  exchanger  can  be  extrapolated  with  Eq.  (24). 

Rdis  =  V(  “  0.0038AV4  +  0.1412AV3  -  1.9399AV2 

+  11. 805AV-  22.023)  (23) 

Rc  =  Rdis4/N  (24) 

4.3.  General  procedure 

The  TSC  is  divided  into  several  design  units  (Du),  each  one 
similar  to  that  presented  in  Fig.  1,  composed  of  two  plates  with  base 
area  of  40  x  150  mm2,  four  thermoelectric  modules  and  one 
dissipater  with  the  corresponding  fan  at  either  side  of  the  plates, 
similar  to  those  described  in  Section  4.2.4.  Therefore,  each  Du 
presents  two  plates  with  circular  channels,  eight  thermoelectric 
modules,  two  dissipaters  and  two  fans.  The  heat  transfer  area  of  a 
Du  is  0.39  m2,  which  corresponds  to  the  external  surface  area  of  two 
dissipaters.  The  TSC  is  composed  of  several  Du,  set  in  series  and 


parallel.  The  design  parameters  and  their  levels  of  variation  are  the 
following  ones: 

Number  of  channels  of  the  HSHE:  12,  36,  60 

Diameter  of  a  channel  of  the  HSHE:  4,  6,  8, 12  mm 

Number  of  Du  set  in  parallel:  1,  2,  3,  4,  5 

The  combination  of  these  parameters  leads  to  60  different  de¬ 
signs  for  the  TSC  to  be  studied.  However,  some  of  them  are  not 
feasible  and  are  not  evaluated.  For  example,  there  is  no  enough 
space  for  60  channels  with  12  mm  of  diameter  in  a  TSC  where  all 
the  Du  are  set  in  series.  Similarly,  a  TSC  with  12  channels  should  not 
have  more  than  2  Du  in  parallel;  otherwise,  the  channels  would  be 
so  separated  from  each  other  that  the  heat  transfer  between  the 
plates  and  the  modules  would  worsen  significantly. 

The  volumetric  flow  rate  is  not  a  fixed  parameter  but  depends 
on  the  design  of  the  TSC.  Therefore,  for  each  of  the  60  combinations, 
four  different  volumetric  flow  rates  are  simulated:  2.22;  2.78;  4.44 
and  6.66  x  10-5  m3/s.  For  each  case,  the  computational  model 
provides  the  number  of  Du  needed  to  evacuate  the  heat  absorbed  by 
the  fluid  in  the  solar  collectors.  Then,  the  real  volumetric  flow  rate 
will  be  the  one  that  meets  Eq.  (5).  Finally,  Eq.  (25)  provides  the 
global  heat  transfer  coefficient  of  the  TSC. 
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U  0.39  Du(Th-Te)  (25) 

5.  Results  and  analysis 

In  order  to  show  the  procedure  followed  to  calculate  the 
working  volumetric  flow  rate  for  each  configuration,  Fig.  4  presents 
that  for  the  TSC  installed  at  the  top  of  the  collectors  ( h  =  0), 
comprising  36  channels  with  diameter  8  mm  and  4  Du  in  parallel. 

If  the  volumetric  flow  rate  in  the  closed  loop  was 
2.78  x  1CT5  m3/s,  Eq.  (7)  would  indicate  that  the  fluid  would  absorb 
2911  W  in  the  collectors  (with  57.9%  of  thermal  efficiency)  and  the 
inlet  temperature  would  be  62.9  °C.  Then,  according  to  Eq.  (6),  the 
thermosyphon  effect  would  provide  a  pressure  increase  of  160  Pa. 
This  point  is  the  second  node  on  the  left  forming  the  solid  curve  in 
Fig.  4. 

Subsequently,  Table  2  would  provide  a  pressure  loss  in  the 
collectors  of  80  Pa  (the  volumetric  flow  rate  in  a  single  collector 
would  be  1.39  x  10-5  m3/s,  since  the  installation  has  two  collectors 
in  parallel).  Likewise,  Eqs.  (9)-(ll )  would  indicate  that  the  pressure 
loss  in  the  connecting  pipes  would  reach  25  Pa,  with  Reynolds  1774, 
friction  factor  0.036  and  16  Pa  of  minor  pressure  losses.  Then,  Eqs. 
(20)-(24)  would  provide  the  thermal  resistances  of  the  heat  ex¬ 
changers  that  the  computational  model  needs  to  indicate  that  the 
TSC  should  include  16  Du  (4  in  series,  4  in  parallel)  to  evacuate  the 
cited  2911  W  of  heat,  and  decrease  the  fluid  temperature  from  90  to 
62.9  °C.  Four  Du  in  series  entail  600  mm-long  channels,  for  a  Rey¬ 
nolds  number  of  107  and  a  pressure  loss  in  the  TSC  of  only  9  Pa, 
according  to  Eqs.  (9) — (11 ).  Therefore,  the  pressure  loss  in  the  closed 
loop  would  reach  114  Pa,  which  is  represented  by  the  second  node 
on  the  left  of  the  broken  curve  in  Fig.  4.  Finally,  Eq.  (25)  would 
indicate  a  global  heat  transfer  coefficient  of  10.43  W/(m2  K). 

This  procedure  is  repeated  for  the  rest  of  the  cited  volumetric 
flow  rates,  and  the  resulting  nodes  form  the  curves  in  Fig.  4.  The 
intersection  of  these  curves  sets  the  working  point,  which  for  this 
configuration  turns  out  to  be  3.35  x  1CT5  m3/s  of  volumetric  flow 
rate  and  130  Pa  of  working  pressure.  All  the  information  regarding 
the  system  performance  for  this  configuration  at  the  working  point 
is  included  in  Table  3,  along  with  the  simulation  results  for  the  rest 
of  the  feasible  configurations.  Those  configurations  with  12  chan¬ 
nels  and  diameter  4  or  6  mm  are  discarded,  since  the  resulting 
pressure  loss  was  so  high  that  the  thermosyphon  effect  does  not 
provide  enough  pressure  to  move  the  fluid.  The  same  occurs  with 
the  configuration  of  36  channels  and  4  mm  of  diameter.  Also,  it  is 


P  (Pa) 


Fig.  4.  Working  point  for  the  TSC  with  36  channels,  diameter  8  mm  and  4  Du  in 
parallel. 


physically  impossible  to  allocate  60  channels  with  diameter  12  mm, 
since  the  cross  area  of  the  channels  is  higher  than  that  of  the  HSHE. 

5  A.  Basic  system  design 

As  can  be  seen  in  Table  3,  the  performance  of  the  whole  system 
follows  logical  patterns.  Firstly,  Fig.  5a  indicates  that  the  total 
pressure  loss  reduces  as  the  cross  section  (S)  of  all  the  channels  of 
the  TSC  increases.  Under  these  circumstances,  the  thermosyphon 
effect  is  able  to  move  higher  volumetric  flow  rates.  Secondly,  as 
Fig.  5b  points  out,  both  the  efficiency  and  the  heat  harvested  by  a 
solar  collector  reduces  as  the  volumetric  flow  rate  increases,  since 
the  average  temperature  in  the  collector  is  higher,  thus  increasing 
the  heat  lost  to  the  environment. 

It  is  common  for  all  thermoelectric  applications,  either  genera¬ 
tors  or  refrigerators,  to  improve  their  performance  when  the 
thermal  resistances  of  the  heat  exchangers  decrease  [16,17  .  In  this 
case,  as  both  thermal  resistances  decrease,  fewer  Du  are  needed  to 
evacuate  the  heat.  Table  3  indicates  that  the  worst  configuration 
requires  20  Du,  whereas  16  are  required  in  the  best  one. 

Regarding  the  HSHE,  since  laminar  flow  must  be  imposed  in  the 
channels,  the  convective  heat  transfer  coefficient  reaches  low  and 
quite  similar  values  in  all  the  cases  [7].  The  thermal  resistance  of 
the  HSHE  decreases  as  the  number  of  channels  increases,  since  the 
total  surface  area  of  the  channels  increases  significantly.  However, 
within  each  value  of  n,  the  thermal  resistance  remains  quite  similar, 
since  the  reduction  in  the  convective  coefficient  is  compensated  by 
the  increase  in  the  surface  area  of  the  channels.  This  fact  indicates 
that  number  of  channels  has  much  more  influence  on  the  thermal 
resistance  of  the  HSHE  than  that  of  the  channel  diameter,  so  it  must 
be  the  main  parameter  to  bear  in  mind  on  designing  the  HSHE  for 
these  applications.  Fig.  6  shows  this  effect. 

Regarding  the  CSHE,  its  thermal  resistance  is  highly  influenced 
by  the  HSHE,  since  the  decrease  in  the  thermal  resistance  of  the 
HSHE  entails  the  increase  in  the  electric  power  provided  by  the 
modules  to  the  fans,  which  rotate  faster,  thus  increasing  the  air 
stream  through  the  fins  and  reducing  the  thermal  resistance  of  each 
dissipater  [8].  However,  this  effect  is  hardly  distinguishable,  since 
forced  convection  is  properly  obtained  in  all  the  cases,  so  the 
thermal  resistance  of  the  CSHE  reaches  similar  values,  as  Table  3 
indicates. 

Finally,  results  report  a  maximum  heat  transfer  coefficient  of 
10.50  W/(m2  K)  for  the  TSC  comprising  60  channels  with  8  mm  of 
diameter,  75%  higher  than  that  provided  by  the  static  dissipate rs 
used  for  overheating  protection  (5.96  W/(m2  I<),  as  Section  2  in¬ 
dicates)  [13].  This  best  configuration  is  composed  of  16  Du ,  4  in 
series  and  4  in  parallel,  with  128  modules  and  8  aluminium  plates, 
each  0.6  x  0.04  x  0.01  mm3  in  size.  The  horizontal  area  occupied 
reaches  0.36  m2,  25%  less  than  that  occupied  by  the  static  dissi- 
paters,  whereas  the  volume  reaches  0.036  m2,  twice  as  high. 

5.2.  System  improvement 

Table  3  indicates  that  the  bottleneck  of  this  application  is  the 
thermal  resistance  of  the  CSHE,  which  is  more  than  three  times  as 
high  as  that  of  the  HSHE  in  the  best  configuration.  In  this  regard, 
Table  4  shows  the  system  performance  for  the  best  configuration  of 
TSC  but  improved  thermal  resistance  of  the  CSHE.  To  conduct  this 
improvement,  Eq.  (24)  is  multiplied  by  decreasing  coefficients,  in 
order  to  reduce  the  thermal  resistance  of  each  dissipater.  Since  only 
the  CSHE  is  modified,  the  heat  absorbed,  collector  efficiency, 
volumetric  flow  rate  and  working  pressure  remains  invariable, 
whereas  the  heat  transfer  coefficient  improves  so  that  fewer  Du  are 
needed  to  evacuate  the  heat. 
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Table  3 

System  performance  for  the  feasible  configurations  of  the  thermoelectric  self-cooling  system. 


n 

D  mm 

S  mm2 

V  x  105  m3/s  QW 

r\  % 

Th°C 

P+bar 

Pc  bar 

Pp  bar 

PjSC  bar  Du  parallel-series 

hh  W/(m2  K) 

Kdis  K/W 

Rh  K/W 

Rc  K/W 

U  W/(m2  K) 

12 

8 

603 

3.08 

2878 

57.3 

76.8 

142 

75 

25 

42 

2  x  10 

289 

0.172 

0.0042 

0.0043 

8.82 

12 

12 

1357 

3.45 

2841 

56.6 

78.5 

125 

86 

30 

9 

2  x  10 

201 

0.173 

0.0041 

0.0043 

8.37 

36 

6 

1018 

3.24 

2858 

56.9 

77.6 

134 

80 

27 

27 

3x6 

322 
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0.0021 

0.0047 

9.56 

36 

8 
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3.35 
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88 

31 

10 

4x4 
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0.0021 

0.0052 

10.43 

36 

12 
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56.4 

78.9 
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89 

31 

2 
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0.0021 

0.0052 

10.34 

60 

4 

754 

2.97 

2912 

58.0 

74.8 

161 

65 

21 

76 

3x6 

283 

0.168 

0.0017 

0.0047 

10.43 

60 

6 

1696 

3.42 

2843 

56.6 

78.4 

126 

85 

30 

11 

4x4 

317 

0.166 

0.0016 

0.0052 

10.50 

60 

8 

3016 

3.51 

2835 

56.5 

78.8 

123 

88 

31 

4 

4x4 

240 

0.166 

0.0016 

0.0052 

10.38 

Q  (w)  n(%) 


Vx  105  (m3/s) 


Fig.  5.  Performance  of  the  TSC  basic  design. 


It  should  be  noted  that  the  original  design  of  the  TSC,  presented 
in  Fig.  3,  includes  a  finned  dissipater  optimized  for  thermoelectric 
applications  [18],  and  it  hardly  exhibits  thermal  resistances  lower 
than  0.15  K/W.  Therefore,  to  obtain  thermal  resistances  significantly 


Fig.  6.  Influence  of  the  diameter  and  number  of  channels  on  the  thermal  resistance  of 
the  hot  side  heat  exchanger. 


lower  than  that,  heat  pipes  must  be  included.  A  heat  pipe  is  a  fin- 
tube  heat  exchanger  that  exhibits  high  heat  transfer  coefficients, 
since  it  includes  a  phase-change  fluid  inside,  leading  to  a  thermal 
resistance  even  an  order  of  magnitude  lower  than  that  of  a  finned 
dissipater  similar  in  size  [19].  A  heat  pipe  presents  improved 
thermal  resistance  when  forced  convection  is  applied  over  the  fins, 
so  the  concept  of  TSC  can  also  be  applied  to  this  heat  exchanger. 

However,  the  last  row  in  Table  4  indicates  that  it  is  useless  to 
reduce  the  thermal  resistance  of  the  heat  pipe  beyond  0.047  K/W, 
since  the  thermal  resistance  of  the  HSHE  becomes  then  the 
bottleneck,  and  the  system  performance  does  not  improve. 
Therefore,  a  TSC  composed  of  12  Du,  with  96  thermoelectric  mod¬ 
ules,  and  heat  pipes  instead  of  finned  dissipaters  would  provide  a 
global  heat  transfer  coefficient  of  13.95  W/(m2  K),  134%  higher  than 
that  provided  by  the  static  dissipaters,  and  33%  higher  than  that  of 
the  previous  design. 

Finally,  as  Section  4.1.1  indicates,  force  convection  allows 
installing  the  TSC  above  the  solar  collector  to  increase  the  pressure 
provided  by  the  thermosyphon  effect,  being  h  the  vertical  distance 
between  them.  Table  5  shows  the  performance  of  the  whole  system 
with  heat  pipes  for  the  TSC,  but  increasing  values  of  h.  Note  that 
36  mm  is  taken  now  as  diameter  of  the  connecting  pipes  to  impose 
laminar  flow  inside  them,  so  the  pressure  loss  in  the  pipes  reduces. 
As  the  vertical  distance  between  the  collectors  and  the  TSC  in¬ 
creases,  the  thermosyphon  effect  provides  higher  volumetric  flow 
rates  at  higher  working  pressure.  It  would  be  represented  in  Fig.  4 
as  if  the  solid  curve  moves  upwards  whereas  the  broken  curve 
remains  constant,  so  that  the  intersection  moves  to  the  right.  As  the 
volumetric  flow  rates  increases,  Reynolds  number  and  the 
convective  coefficient  increases,  according  to  Eqs.  (11)  and  (20), 
leading  to  a  final  increase  in  the  global  heat  transfer  coefficient. 
However,  given  that  laminar  flow  is  still  imposed  in  the  channels  of 
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Table  4 

System  performance  for  the  best  configuration  of  the  thermoelectric  self-cooling  system  with  improved  thermal  resistance  in  the  cold  side  heat  exchanger. 


V  X  105  m3/s 

Q  W 

H% 

Th°C 

P>  bar 

bar 

Pp  bar 

PTSc  bar 

Du  parallel-series 

hh  W/(m2  K) 

Kdis  K/W 

Rh  K/W 

Rc  K/W 

U  W/(m2  K) 

3.42 

2843 

56.6 

78.4 

126 

85 

30 

11 

4x4 

317 

0.166 

0.0016 

0.0052 

10.50 

3.45 

2841 

56.6 

78.5 

126 

86 

30 

10 

5x3 

350 

0.099 

0.0016 

0.0033 

11.16 

3.45 

2841 

56.6 

78.5 

125 

86 

30 

9 

4x3 

350 

0.047 

0.0020 

0.0020 

13.95 

3.45 

2841 

56.6 

78.5 

125 

86 

30 

9 

4x3 

350 

0.023 

0.0020 

0.0009 

13.95 

Table  5 

System  performance  for  the  best  configuration  of  the  thermoelectric  self-cooling  system  with  improved  thermal  resistance  in  the  cold  side  heat  exchanger  and  increased 
thermosyphon  effect. 


h  m 

V  x  105  m3/s 

Q  W 

V% 

Th°C 

P1  bar 

Pc  bar 

Pp  bar 

PjSC  bar 

Du  parallel-series 

hh  W/(m2  K) 

Kdis  K/W 

Rh  K/W 

Rc  K/W 

U  W/(m2  K) 

0 

3.45 

2841 

56.6 

78.5 

125 

86 

30 

9 

4x3 

350 

0.023 

0.0020 

0.0009 

13.95 

0.5 

4.56 

2774 

55.3 

81.7 

146 

125 

14 

7 

5x2 

444 

0.022 

0.0021 

0.0011 

15.23 

1 

5.18 

2748 

54.7 

82.9 

174 

148 

18 

8 

5x2 

464 

0.022 

0.0020 

0.0011 

14.73 

the  TSC,  this  effect  has  limited  impact  on  the  system  performance, 
and,  for  h  higher  than  0.5,  it  becomes  so  insignificant  that  no 
additional  DU  can  be  removed. 

Therefore,  the  best  design  of  a  TSC  to  provide  overheating  pro¬ 
tection  to  the  real  solar  installation  of  this  work  would  be  elevated 
0.5  m  over  the  top  of  the  collectors,  and  would  be  composed  of  10 
Du,  with  80  thermoelectric  modules,  and  heat  pipes,  for  a  global 
heat  transfer  coefficient  of  15.23  W/(m2  K),  155%  higher  than  that 
provided  by  the  static  dissipaters  and  9%  higher  than  that  of  the 
previous  design. 

6.  Conclusions 

The  main  conclusion  of  this  paper  is  that  it  is  feasible  to  develop 
a  zero-power-consumption  system  to  protect  solar  collectors  from 
overheating.  The  combination  of  thermoelectric-self-cooling  tech¬ 
nology  and  thermosyphon  effect  provides  a  device  that  dissipates 
the  excess  heat  and  requires  no  electric  power,  thus  being  able  to 
work  properly  under  the  most  unfavourable  conditions,  such  as 
electric  power  failures,  circulator  pump  failures,  maintenance  shut¬ 
downs,  etc. 

The  simplest  design  of  the  thermoelectric-self-cooling  system 
includes  128  thermoelectric  modules  Kryotherm  TGM-287-1.0- 
1.5,  and  the  corresponding  heat  exchangers.  The  hot  side 
heat  exchanger  is  composed  of  8  aluminium  plates, 
0.6  x  0.04  x  0.01  mm3  each,  including  60  circular  channels  with 
8  mm  of  diameter,  through  which  the  fluid  circulates.  The  cold  side 
heat  exchanger  is  composed  of  32  finned  dissipaters,  each  one 
comprising  a  square  base  plate  with  side  length  150  mm  and 
height  12  mm,  and  23  fins  with  dimensions  150  x  25  x  1.5  mm3. 
Finally,  32  fans  Sunon  KDE1208PTS1,  energized  by  the  modules, 
provide  forced  convection  over  the  dissipaters. 

When  the  fluid  temperature  at  the  exit  of  the  collectors  reaches 
90  °C,  the  fluid  is  forced  to  flow  through  the  closed  loop  composed 
of  the  collectors,  the  thermoelectric-self-cooling  system  and  the 
connecting  pipes,  where  the  thermosyphon  effect  creates  the  force 
required  to  make  the  fluid  circulate.  The  hot  fluid  in  the  channels  of 
the  hot  side  heat  exchanger  induces  a  temperature  difference  be¬ 
tween  the  ends  of  the  thermoelectric  modules,  which  provide 
electric  power  to  the  fans.  Then,  forced  convection  is  obtained  over 
the  dissipaters,  leading  to  a  heat  transfer  coefficient  of  10.50  W / 
(m2  I<)  without  electricity  consumption,  75%  higher  than  that 
provided  by  the  static  dissipaters  currently  used  for  overheating 
protection. 

The  thermal  resistances  of  the  heat  exchangers  emerge  as  key 
parameters  in  the  system  performance,  so  their  improvement  in¬ 
creases  the  global  heat  transfer  coefficient.  Thus,  the  inclusion  of 


heat  pipes  instead  of  dissipaters  in  the  cold  side  heat  exchanger 
leads  to  a  heat  transfer  coefficient  of  13.95  W/(m2  I<),  33%  higher 
than  that  of  the  first  design,  which  allows  a  25%  decrease  in  the 
number  of  modules,  from  128  to  96.  Finally,  the  improvement  of  the 
thermal  resistance  of  the  hot  side  heat  exchanger  comes  from  the 
increase  in  the  vertical  distance  between  the  dissipation  system 
and  the  top  of  the  collectors,  which  boosts  the  thermosyphon  ef¬ 
fect,  allowing  higher  volumetric  flow  rates  at  higher  working 
pressures  and,  as  a  consequence,  higher  convective  heat  transfer 
coefficients  between  the  fluid  and  the  channels.  This  measure  en¬ 
tails  an  additional  9%  increase  in  the  global  heat  transfer  coefficient, 
reaching  15.23  W/(m2  I<),  and  a  17%  decrease  in  the  number  of 
modules  required,  from  96  to  80.  Compared  to  static  dissipaters,  the 
global  heat  transfer  coefficient  is  improved  by  155%. 

To  conclude,  the  system  joins  the  advantages  of  dynamic  and 
static  dissipaters,  since  it  overcomes  the  restrictions  related  to 
natural  convection  and  consumes  no  electricity.  Also,  it  presents 
the  robustness  and  reliability  exhibited  by  all  thermoelectric  ap¬ 
plications,  thus  standing  out  as  a  promising  alternative  to  prevent 
the  overheating  of  solar  collectors. 
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